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ABSTRACT 



Measurements of black hole spin made using the continuum-fitting method 
rely on the assumption that the inclination of the black hole's spin axis to our line 
of sight is the same as the orbital inclination angle i of the host binary system. 
The X-ray and radio jet data available for the microquasar XTE J1550-564 offer 
a rare opportunity to test this assumption. Following the work of others, we 
have modeled these data and thereby determined the inclination angle 9 of the 
jet axis, which is presumed to be aligned with the black hole's spin axis. We find 
9 ~ 71° and place an upper limit on the difference between the spin and orbital 
inclinations of \9 — i\ < 12 deg (90% confidence). Our measurement tests for 
misalignment along the line of sight while providing no constraint perpendicular 
to this plane. Our constraint on the misalignment angle supports the prediction 
that the spinning black hole in XTE J1550-564 has aligned itself with the orbital 
plane and provides support for the measurement of its spin via the continuum- 
fitting method. Our conclusions are based on a simple and reasonable model of a 
pair of symmetric jets propagating into a low density cavity whose western wall 
is ~ 20% closer to XTE J1550-564 than its eastern wall. 

Subject headings: black hole physics — stars: individual (XTE J1550-564) - 
X-rays: binaries 



1. Introduction 



Although it is thought that the Galaxy is host to tens of millions of stellar-mass black 
holes, only about 50 have been discovered (jOzel et al.l |2010| ). All of them are accretion- 
powered X-ray sources that are located in X-ray binary systems. Most such systems, which 
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are similar to our featured black hole binary XTE J1550-564, have short orbital periods 
(P ~ 1 d) and are comprised of a low-mass (< 1 M ) donor star and a ~ 10 M black hole. 
A stream of gas from the Roche-lobe-filling star feeds into the outer part of an accretion disk 
that encircles the black hole. On a time scale of weeks, viscous forces in the disk cause this 
gas to move radially inward to the center. Within a few hundred kilometers of the black hole, 
the optically-thick gas reaches a temperature of ~ 10 7 K and produces an X-ray luminosity 
that is near the Eddington limit (L ~ 10 39 erg s _1 ). Accretion onto the black hole is not a 
steady process: A typical source is luminous for only about a year, and then it fades into a 
quiescent state for years or decades. 

XTE J1550-564 (hereafter J1550) is a much-studied Galactic black-hole transient sys- 
tem that was discovered on 1998 September 6 using the All-Sky Monitor (ASM) onboard the 
Rossi X-ray Timing Explorer (RXTE). Thereafter, it was obser ved almost daily dur ing its 
entire 8- month outburst cycle using RXTE's pointed instruments ( jSobczak et al.ll2000l ). Two 
weeks into outburst, the source abruptly rose fourfold in intensity and produced a brilliant 
7-Crab fl are. During this X -ray flare, J1550 was approximately at its Eddington limit for 
~ 1 day ( jSteiner et al.ll201ll ). Four days later, radio observations made using the Australian 
Long Baseli ne Array (LBA) revealed relativistic ejecta moving both eastward and westward 
from J1550 (IHannikainen et al.ll2009l ). The two components were observed to be separated 
by ~ 250 mas and moving at relative speed of /i app ~ 65 mas/d, equivalent to an apparent 
separation velocity of ~ 1.7c. Nearly two years later, Chandr a imaging observat ions revealed 
large-scale (> 20") relativistic jets undergoing deceleration (jCorbel et al.l 120021 ). This land- 
mark discovery of a pair of ballistic X-ray jets was the first detection of its kind for a Galactic 



source. 



By modeling an extensive collection of optical and infrared data for J1550. lOrosz et al. 

(120111 ) have determined the mass of the black hole primary, M = 9.1±0.6 M , the distance to 
the binary, D = 4.38^Q 4^ kpc, and the inclination of its orbital plane, i = 74?7±3?8. Assum- 



i ng th at the black hole's spin is aligned with the orbital angular momentum, ISteiner et al. 
( 1201 ll ) have measured the spin using the continuum-fitting method to be a* = 0.34 



+0.20 
0.28' 



where a* = cJ sp - m /GM is the black hole's dimensionless spin parameter and J sp i n its angular 
momentum. Steiner et al. also measured the spin using the independent Fe-line method and 
find a* = 0.55^1^; taken together, the two measurements imply a* « 0.5. The continuum- 



fitting method relies on a model for the thermal emission from an accretion disk (jZhang et al. 



19971 ). while the Fe-line method rel ies on a model of the relativistically broadened fluores- 
cence features emitted by the disk (IFabian et al.lll989l ). 
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For a black-hole binary system like J1550, with a low-mass companion, the ratio of the 
orbital angular momentum to the spin angular momentum of the black hole is given by 

JorW J spiD - 65 a; 1 ( JL-) , (1) 



v 10 M Q J \M q J \1 dj 

where M2 is the mass of the secondary star. For J1550, this ratio is ~50, and thus it 
is reasonable to expect that, given a means of interaction, the spin of the black hole will 
eventually come into alignment with the orbital angular momentum. The time scale for 
this to occur is an important question for continuum-fitting spin measurements because in 
applying this method one generally must assume that the two vectors are aligned. 

If there is an initial misalignment between the spin and the orbital angular momenta, 
then Lense-Thirring precession will cause the inner X-ray-emitt ing portion of the disk to line 
up with the spin of the black hole ( IBardeen &: Pettersonlll975l ). At the same time, at very 
large scales, the disk will align itself with the orbital plane, and the transition between these 
regimes will manifest as a warp in the disk. When a misalignment is present, the black hole 
will be torqued in to alignment by the accreti ng matter acting with a lever arm of order the 



warp radius (e.g., iNataraian fc Pri ngle 



1998). Using a maximally conservative (minimum- 



torque) assumption, iFragos et al.l ( 1201 01 ) concluded (based on a population synthesis study) 
that the spin axes of most black hole primaries will be tilted less than 10°. Fragos et al. 
assumed that the torque acts at the innermost stable circular orbit, -Risco < 6GM/c 2 for 
a* > 0, whereas the warp radius has be en estimated to be located at i? w m 200GM/c 2 
king et al-lliooi lLodato fc Pringlelhood ). 



For a typical system, the time scale for accretion to torque the black hole into alignme nt 



has been estimated to be £ a iign ~ 10 6 — 10 s years ( Martin et al. 2008 ; Maccarone 2002 )^1. 



Therefore, one expects alignment to occur early in the lifetime of an old-population transient 
system, such as J1550, and that most such systems will presently be well aligned. 

It is obviously important to test this theoretical expectation. However, it has proved 
challenging to obtain a firm measurement of the degree of alignment for any black hole binary. 
Such a measurement requires a determination of the position angle of the binary on the plane 
of the sky. While this may be possible in the future for J1550, we lack the requisite orbital 
astrometric data for the system and are therefore limited to testing for alignment along the 
line of sight. Measuring the orbital inclination angle of the binar y is relatively simple and is 



routinely done by modeling optical data (e.g.. lOrosz et al.l 120091 ). In contrast, it has proved 



difficult to obtain reliable estimates of the inclination of the inner disk. 



1 Maccarond ( 2002 ) overestimated t a i; gn as the result of a numerical error in his Eqn. 6, which implies a 
time scale that is 50 times longer than that implied by his Eqn. 1. 
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Currently, the most direct way of determining the inner-disk inclination is by modeling 
jet ejecta, which are presumed to be ali gned with the black hole's s pin axis. For the case 
of symmetric ejecta, see the review by iMirabel &: Rodriguezl (119991 ) . The jet ejecta that 
are relevant to this paper are pairs of discrete, detectable condensations of radio-emitting 
plasma, which we generally refer to in shorthand as "jets." 

An alternative approach to measuring the inner-disk inclination is y ia the same Fe-line 
method that is used to measure black hole spin ([Reynolds fc Nowakll2003l ). However, existing 
models make simplifying assumptions concerning how the ionization state of the disk varies 
with radius. Given that there is a degeneracy between ionization and inclination in Fe- 
line/reflection models for stellar- mass black holes, these inclination estimates are subject to 
a systematic uncertainty of unknown magnitude. Meanwhile, prospects are good that more 
advanced reflection models will provide robust estimates of inclination. 

Based on observations of radio jets, two confirmed black hole systems, GRO J1655-40 
and SAX 1819-2525, are good candidates for hosting misaligned black holes. In the case of 
GRO J1655-40, using a k inematic model for the jets and measurements of proper motion, 
Hjellming fc Rupenl (119951 ) reported a jet inclination angle of 85°. However, the authors give 
no error estimate for either the jet inclination angle or the proper motion. Furthermore, the 
reliability of the estimate for the jet inclination angle is called into qu estion by the intrinsic 
and variable asymmetries that were observed for the opposing jets (IMirabel fc Rodriguez 
19991 ). Taking the 85° jet inclination an gle at face value, on e concludes that the jet axis and 
orbital vector are misaligned by > 15° (IGreene et al.ll200ll ). 



In the case of SAX J1819-2525, the evidence is less certain. There is only a single 
observation of extended radio emission (because the source faded promptly). By making the 
assumption that this emission was associated with a major X-ray outburst that occurred 
hours ea rlier, superluminal motion (0 app > 10c) and a misalign ment angle of > 50 ° were 
inferred (jOrosz et al.l l2001t iHjellming et al.ll2000f ). However, as IChaty et al.l (120031 ) have 
argued, the jet may have been ejected a couple of weeks before the major outburst, in 
which case the Lorentz factor of the jet was modest and its inclination consistent with the 
inclination of the binary. This is a reasonable possibility given that the source was observed 
to be active at optical wavelengths for several weeks before the X-ray outburst. 

Compared to the jets in these two systems and those in other Galactic microquasars, the 
jet ejections observed for J1550 are remarkable. They were observable for years (rather than 
weeks or months), and therefore their physical separation from J1550 was observed to become 
exceptionally large. These are possibly the largest resolved jets observe d for any black hole 



when c onsidering; th e dimensionless distance between them, i.e., d/M ( iHao fc Zhang! 12009 



and see iHeinzl |2002| ). By this measure, the maximum 0.7 pc distance between the jet and 
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J1550 corresponds to 7 Mpc for a supermassive black hole of 10 8 M Q 



In a previous study, IWang et al.l ( 120031 ) modeled the evolution and light curve of J1550's 
ballistic jets using the same model we employ, namely, an expanding jet interacting with the 
interstellar medium (ISM). They modeled the data for the eastern jet, attributing the X-ray 
emission to a reverse shock, and fo und that the gas dens ity around J1550 is unusually low. 
Later, their work was extended by lHao fc Zhangl (120091 ) to include the western jet. Both 
groups focused their attention on the properties of the environment around J1550; accord- 
ingly, they adopted nominal and fixed values for jet inclination (50° and 68°, respectively), 
initial Lorentz factor (3), and jet energy (3.6 xlO 44 erg). Both groups found evidence for the 
existence of a low density cavity around J1550 (modeled in more detail by Hao & Zhang), 
and a possible east-west asymmetry in the ambient gas. 

While we follow in the footsteps of Wang et al. and Hao & Zhang, our aim is different. 
We are focused on the question of the alignment of the inclination angle of the black hole's 
spin axis and the orbital inclination angle. Therefore, in distinction with the earlier work, we 
disregard the X-ray light-curve data, which are primarily useful in constraining the emission 
mechanisms or the electron density and magnetic fields in the jet. Rather, we concentrate on 
modeling the kinematics of the ballistic jets and deriving reliable values and error estimates 
for the kinematic parameters. The parameter of chief interest is the inclination of the black 
hole's spin axis. 



2. Data 

We use archival Chandra X-ray Observatory data for eight observations of J1550 that 
were obtained using the Advanced CCD Imaging Spectrometer (ACIS) between 2000 June 
and 2003 October. The exposure times range from 4 — 50 ks. Pipeline processed level-2 
event filea^] were used to produce images of the field of J1550. When detected, images of 
the eastern (approaching) jet yielded 16-40 counts and the western (receding) jet 100-400 
counts; J1550 itself was always detected and yielded 60-3000 counts. 



These same Chandra data were used by lHao fc Zhangl (120091 ) in their analysis of the 



X-ray jets. They relied on the absolute astrometric precision of Chandra in order to derive 
positions for each jet and thereby its offset from J1550. We have reduced the astrometric 
errors severalfold by directly measuring in each image the relative separations between J1550 
and the jets. 



2 using CXC DS-7.6.10 



-6 - 



In measuring the precise jet positions, which are given in Table [U we smoothed each 
image using a 1" Gauss ian kernel and then determined the centroid of each jet using the 
DAOphot find routine (IStetsonlll987l ). This procedure was used to derive initial estimates 
for all the jet positions. Then, 1000 Poisson random realizations of each field were produced, 
and the centroid measurements were repeated. In most cases, the positions for a given jet 
were tightly clustered about a single value, and a separation and error were derived from 
this distribution. However, for three observations of the eastern jet (Obs. XI, X3, and 
X6 in Tabled]) the images are particularly faint (possibly because the emission is extended), 
which resulted in a broad distribution of positions. In these Gaussian-weighted mean 

based on the jet position angle (f)j for each realization j was used to derive the separation 
between the jet and J1550. As a reference value, we used the average position angle for the 
jets 4>pa along with its error opA, 4>pa = 94?25 ± 0?3 (mea s ured e ast o f north). Thi s value 
is consistent with those determined by lHannikainen et al.l (120091 ) and ICorbel et al.l (120021 ) 
and was measured for a single frame generated by coaligning and coadding all of the X- 
ray images. The weights Wj were calculated according to \og(wj) = —\{4>j — 0pa) 2 /°pa- 
Typically, the position errors for the eastern jet were several tenths of an arcsec, while for 
the brighter western jet they were < 0.1". 

In addition to the positions derived using the Chandra data, we include in our analysis 
two radio positions (Obs. Rl and R2 in Table [1]). These measurements are taken from 
Corbel et al.l (120021 ). who derived positions from observations obtained using the Australia 
Telescope Compact Array (ATCA) on 2000 June 1 and 2002 January 29. In the first obser- 
vation, only the eastern jet is observed, whereas in the second, the eastern jet has faded and 
the western jet alone is present. 

As a final constraint on our kinematic jet model, we require that the apparent sepa- 
ration speed of the jets a t laun ch match the value measured using the LBA, 65.5 ± 13.2 
mas/d lHannikainen et al.l (120091 ) . This speed and the jet positions are the sole inputs to our 
principal model in Section [51 However, in Sec tion |6l we additionally consider radio inten- 
sity measurements given by lHannikainen et al.l (120091 ) . They report 2.29 GHz flux densities 
taken four and six days after the X-ray flare with intensity ratios of Sei/Swi = 3.55 and 
Se2/Sw2 = 2.40, respectively; we assume that these ratios are uncertain by 25%. We also 
adopt their measurements of the radio spectral index, ipi = —0.43 and ip2 — —0.21, taken 
from flux densities measured with the ACTA at 4.8 and 8.6 GHz. The spectral index mea- 
surements and corresponding LBA images, while not strictly simultaneous, were obtained 
within several hours of one another. 
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The Jet Model 



The development of our kinematic jet model follows lHao fc Zhana ( 120091 ) and lWang et al. 



( 120031 ). The model we use has been designed to describe gamma-ray bursts, but it is applica- 
ble to a relativistic, adiabatically expanding jet. To begin, we consider a pair of symmetric 
jets, each launched with a kinetic energy E and Lorentz factor r . As the jets expand 
into their environments, they entrain material from the surrounding medium, dissipate their 
kinetic energy at the shock front and heat the ISM. We negl ect radiative losses and assume 
that the jets are confined and evolve adiabatically. Following [Wang et al.l ( 120031 ). we assume 
that particles are accelerated uniformly and randomly at the shock front. Each such jet 
obeys the relation 

Eo = (T- l)M c 2 + a{T 2 sh - l)m sw c 2 , (2) 

where T is the instantaneous bulk Lorentz factor of the jet, Mq is the mass of the jet ejecta, 
and T s h is the Lorentz factor at the shock front. The mass of the entrained material, m sw , that 
has been swept up by the shock is approximately m sw = 6 2 m p n7r_R 3 /3, where 6 and R are 
respectively the jet half opening angle and the distance the jet has traveled. The numerical 
factor a varies from Q.35 for ultrarelativisti c shocks to ~ 0.73 in the nonrel ativistic limit 
jWang et all 1200.4 blandford fc McKee! Il976l ). Following Euang et all Jl999h . we adopt a 
simple numerical scaling to interpolate between the two regimes: a = 0.73 — 0.38/3, (/3 = 

At the shock fr ont, the jump condition re lates the bulk Lorentz factor of the jet to that 
of the shocked gas (jBlandford fc McKedll976l ): 



F 2 

1 sh 



(r + i)[7(r -i) + i]- 

7(2-7)(r-l) + 2 



(3) 



The adiabatic index 7 varies between 4/3 and 5/3, which are respectively its ultrarelativis- 
tic and nonrelativi s tic limits. We interpolate between th ese regimes via 7 = (4r + l)/3r 
jrluang et al"1ll999[ IWang et aI.ll20o4 lrIao fc Zhan~g1l2009h . 



On the plane of the sky, the apparent proper motions of the approaching and receding 
jets, \i a and /i r , are given by 

(3 c SU16 1 (3 c sin8 



Ha 



D(l - (3 cost 



D(l + j3 cos6 



(4) 



As we show in Section [51 the simple mo del g overned by E qn. [2]fails to fit the observations. 
Motivated by the results of lHao fc Zhangl and IWang et al.l . we have generalized Eqn. [2] to 
allow for the jets to first propagate through a low density cavity before encountering and 
shocking against the ISM. In the east-west direction, we allow for the cavity to differ in 
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size. We additionally consider the possibility of an intrinsic asymmetry in the jets. Eqn. [2] 
becomes: 

V E = (r - l) V M c 2 + a(T 2 h - l)m sw c 2 , (5) 
and the entrained mass is now 



Q 2 m p nir [ R 3 , R< (R 



° X ^ (CR C r) 3 + 5[R 3 ~ (CRcr) 3 ], R>(rZ (0) 



where R cr and 5 are respectively the radius of the cavity centered on J1550 and the density 
jump at the cavity boundary. The ratio of the western-to-eastern cavity dimensions is given 
by (. Similarly, 77 = (£°/^02) WCS t t parameterizes the asymmetry of the jets. In application, the 
asymmetry parameters ( and rj are taken to be unity for the eastern jet and can vary for the 
western jet. 

In order to obtain a model solution for a particular set of parameters, we evolve the 
energy equation as the jet expands (either Eqn. [2] or Eqn. [5]) in 4-hour time-steps by se- 
quentially solving for T(t) in the rest frame of J1550. At each time step, we calculate the 
separation between each jet and the central source by integrating f3(t) and by calculating 
the projected angles a: a(t') = R(t)sm 9/D. Here, t' = t =F R(t) cos 9/c is the observer's 
time, which takes into account the time delay between J1550's rest frame and that of the 
observer for whom the light-travel paths of the approaching and receding jets are respectively 
shortened and elongated. 

Our model requires up to eight physical parameters: 9, T Q , D, R cr , 5, r], (, and lastly 
the "effective energy" E which we now define. As alluded to above, a degeneracy exists in 
our model between jet energy, ambient gas density, and the jet opening angle. These three 
quantities appear as a single and inseparable term in the kinematic equations, Eq/hQ 2 . To 
make physical sense of this combined quantity, we assume that the de nsity of the ISM is a 



standard nisM = 1 cm 3 , so that n = 1/5 cm 3 , and adopt 6 = 1° ( jKaaret et al.ll2003l ). 



Predicated upon our assumed values for n^M and 0, the jet energy E is then E = E. 

Finally, we go beyond our principal, kinematic model to consider the ratio of the radio 
intensities of the two jets. We consider the simplest case of the ejection of a pair of identical 
and unimpeded condensations. When measured at equal separation from the black hole, one 
has 

S a fl + (3cos9^ 3 ^ 



S r \1- (3 cos 9 J ' ^ 
where ip is the spectral index and the subscripts a and r refer to the approa ching and 



receding jets, which are taken to be discrete ejecta (IMirabel fc Rodriguez! 1 1 99 9l ). Because 



the jets are observed at unequal distances from the black hole, we must adopt a model of how 
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jet intensity varies with time; we assume a simple power-law dependence. Then, allowing 
for our case of intrinsically asymmetric jets, Eqn. [7] becomes 

Sa 

S r 

where A is a fit parameter, which for positive values describes a decay in brightness with 
time. The effect of jet asymmetry on the radio emission is captured by q, which can range 
from -1 to 1.5 depending on the source of asymmetry: E (q = —1), n (q = 0), or © 
(qe [1,1.5]). 



T r (l + j3 r cos 

r a (l - P a COS 



O— : tp — £A 



4. Markov Chain Monte Carlo 

Markov Chain Monte Carlo (MCMC) is a powerful statistical technique by which ran- 
dom samples are drawn from a posterior distribution of arbitrary form. In our case, the 
posterior distribution is the probability of our model parameters, given the data. MCMC 
algorithms perform a "guided walk" of transitions through parameter space such that, after 
an initial burn-in phase, the chain directly reproduces the likelihood surface for the model. 
MCMC has several advantages over traditional gridded-search algorithms when the num- 
ber of parameters is large. For example, the search time with MCM C scales approximatel y 



linearly with the number of parameters rather than exponentially ( iMartmez et al.ll2009l ). 
Furthermore, the ergodic property of the Markov Chain guarantees (asymptotically) that 
the chain will fully explore parameter space and reach the optimum global solution. 

Transitions in the chain are effected via a "jump" distribution^ J(x*\x n ) (e.g., a mul- 
tivariate Gaussian) that defines a probability of selecting a candidate transition to a new 
state x* given the current state x n . T he transition probability from x n to x* is governed 



by the Metropolis-Hastings algorithm ([Hastings! 119701 ) and is determined by the ratio r of 
probability densities 

p(x*\y)J(x n \x*) 
p(x n \y)J(x*\x n y 

where y refers to the data, and p(a\b) should be read in the usual way as the probability of 
a given b. The term p(x*\y)/p(x n \y) in the equation above gives the likelihood ratio of the 
two states, while the remaining term corrects for bias introduced by the jump-distribution 



3 We implement a particular class of the algorithm known as random- walk MCMC. In this approach, a 
sequence of transitions from the current parameter values are proposed and are then incrementally accepted 
or rejected. 
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density at each state. The state of the next link in the chain, x n +i, is then chosen according 
to 

f x*, with probability min[r, 1], 
\ x n , otherwise. 

The likelihood ratio appearing in Eqn. [9] is calculated by evaluating the \ 2 for each state 
while taking into account the prior p on all of the model parameters. In this case, the priors 

N 

are introduced independently so that p = Y\ pk, where N is the number of parameters and 

fc=i 

pk gives the prior for parameter k. Omitting additive constants, the log-likelihood for state 
x is 

log {p{x\y)) = -- [ X 2 (x) - 2 log {p{x))\ . (11) 



4.1. MCMC in Practice 

As Eqn. [TT] makes apparent, the prior acts as a penalty to x 2 , and for the special case 
that a prior is "flat" (i.e., independent of x), one recovers the usual least-squares formula. 
It is also worth noting that because the prior only enters into the MCMC chain generation 
as a ratio (Eqn. [9]) the scaling of the prior is arbitrary. We introduce a new term for this 
penalized x 2 , namely x 2 , sucn that 

xl(x) = X 2 (x)-2log(p(x)). (12) 

Unless stated otherwise, we choose to normalize the prior so that the penalty term (2 log[p(x)]) 
is zero at the best fit, i.e., at the minimum value of x%- 

We adopt an asymmetric Gaussian prior on the distance becaus e it has been previ ously 



measured using optical and near-infrared photometry (see Sectional lOrosz et al.ll201ll ). For 
the asymmetric terms, we adopt a log-flat prior on the difference from unity, i.e., p v oc 
min[l / 77, 77] (and likewise for p^). As an example and stated differently, we consider a term 
implying a 10-fold asymmetry to be a priori one tenth as likely as one that is symmetric. We 
adopt flat priors on 9 and R CI and flat priors on the log- values of scale parameters (i.e., the 
jet energy, T and S). The priors and parameter rangesj are discussed further and illustrated 
in Section O 



4 While it is optimal to use an unbounded parameter space in performing MCMC sampling, it is also 
sensible to set physically meaningful constraints on the parameters (e.g., Tq > 1). To achieve both objectives, 
we have transformed each parameter using a logit function to map a parameter z from its range [z m i n , z m ax] 
onto an infinite scale: logit(t) = z mnl + (z max — z m - m )/(l + e~*) for — 00 < t < 00. 
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In order to initialize the chain and the jump distribution, we make starting guesses for 
the model parameters and step sizes. These initial values are improved upon by running 
a sequence of "training" iterations. The training phase incrementally improves the jump 
function until its shape is a close approximation to the posterior covariance matrix, thereby 
greatly increasing the MCMC efficiency The sequence becomes increasingly tuned to the 
likelihood surface, simultaneously refining E (the covariance estimatejj and optimizing the 
solution. 

The training phase continued for a minimum of 15 iterations, each of which generated 
a trial chain with 2000 elements. Training terminated either after 25 cycles were completed 
or when the chain attained an acceptance fraction between 24% and 37%B 

Upon completing the training cycle, 8 chains were generated and run in parallel using 
the trained jump function, each to a length of 110 thousand elements. Seven of the starting 
positions were chosen by sampling using a dispersed covariance E' = 10 E about the final 
training position, and the eighth was started directly from the end location reached by the 
training sequence. The initial 10,000 elements of each chain were rejected as the "burn-in" 
phase during which the chains relax toward a stationary distribution. Our final results are 
based on a total of 8 xl0 5 MCMC samples. Convergence of the MCMC run is determined 



using the criterion of iGelman fc Rubin! (119921 ). R. The closeness of this criterion to unity is 



the measure of convergence. 

In Figure dj we plot a trace of our parallel runs over time for inclination in our adopted 
model (see Section [5]). In the bottom panel, we show the Gelman & Rubin convergence 
diagnostic of the chain oyer tim e. Typically, a chain is considered converged if R < 1.1, or 
1.2 (see, e.g., Verde et al. 2003 For 9, our parameter of interest, we obtain R < 1.01. 



£ is calculated from the chain positions and is used to define the jump function for each sequence. The 
jump function is taken to be a t— distribution with 4 degrees of freedom that is symmetric about the present 
position. 

6 The target acceptance fraction was set at ss 32%. The o ptimal value ranges f rom 23% for an infinite- 
dimensional problem to « 45% for a univariate problem ([German et al.lll996T ). Each run produced an 
acceptance fraction of at least 20% 

7 Larger values of R suggest that either the parameter space is insufficiently sampled or that the chains 
are not fully evolved. 
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Table 1. Relative Jet Positions 



Obs. 




Eastern Offset (arcsec) 


Western Offset (arcsec) 


Rl 


620.5 


21.9±0.3 b 




XI 


628.5 


21.5±0.5 




X2 


701.4 


22.7±0.2 




X3 


722.2 


23.7±0.5 




R2 


1227.5 




22.6±0.3 b 


X4 


1268.8 


28.5±0.2 


22.78±0.05 


X5 


1368.5 




23.19±0.07 


X6 


1466.0 


29.6±0.6 


23.44±0.10 


X7 


1591.3 




23.76±0.10 


X8 


1859.6 




24.4 ±0.2 



a Time since the jets were launched: MJD - 51076. 



b ACTA position from lCorbel et all (120021 ). 




Fig. 1. — top: The trace of 9 for Model AC of Section [51 Eight parallel chains are used; 
for each, the initial 10 4 elements are generated during the burn-in phase and discarded from 
the analysis, bottom: The convergence of the chain over time. The chains reach convergence 
quickly, which is indicative of efficient sampling. 
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5. Results 

In this section we consider three symmetric-jet models, including our adopted model. 
For these models, and for the additional models discussed in the following section, we assume 
that the jets were launched at the time of J1550's giant X-ray flare (Section [TJ. 

5.1. Two Preliminary Models 

We first consider and rule out two simple models. For the simpler of these, which we 
refer to as Model SI, the jets are symmetric and propagate through a uniform medium 
(Eqn. El i.e., r) = ( = 5 = 1 and R cr = 0). The strong deceleration of the jets at late times 
is not accommodated by this model, and the best fit achieved is unacceptable, X 2 p jv — 68. 
For Model S2, we introduce a symmetric cavity centered on J1550 with S and R CT as free fit 
parameters. The fit is significantly improved, y^ p j v = 42, but it is still far from acceptable. 
The results for both models are given in Table [21 

5.2. Our Adopted Model 

We now consider our primary model - an extension of Model S2 that allows the source to 
be positioned off-center in the cavity. This asymmetric cavity model (Model AC) is obtained 
by freeing the fit parameter ( (while leaving 7] fixed at unity; see Eqn. [5]). As illustrated in 
Figure El for a modest (22%) degree of asymmetry, this model produces a successful fit to 
the data with x%l v = 1-44. Results are given in Tableland marginal distributions from 
the MCMC run are shown for each parameter in Figure [31 

The eastern and western cavity walls are located respectively at 0.6 pc and 0.5 pc from 
the black hole and the density contrast at the boundary of the cavity is ~ 100. The gas 
density within the cavity is much lower than that of the ISM. This must be the case in order 
for the jets to have passed through without sweeping up enough mass to halt their expansion. 
Motion within the cavity lasted for f=al.5 years (in the frame of J1550), until the receding 
western jet impacted the dense ISM at its cavity wall and abruptly began decelerating in 
advance of its eastern counterpart (see Fig. [2]). 

The total energy for both jets is an impressive E tot 10 46 erg -^-3 ( i deg ) " ^ l auncn ; 
the Lorentz factor of the jets is constrained to be T > 1.6 (99.7% confidence). However, the 
data provide no upper limit on T, as implied by Figure El which shows that for large values of 
T the distribution closely tracks the prior. Likewise, the data only weakly constrain J1550's 
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distance. However, the remaining five parameters are well determined by the data and are 
quite independent of their priors (Fig. [3]). For the key parameter, the jet inclination angle, 
we obtain 9 ~ 71° (64° < 9 < 83° at 90% confidence) and find only moderate correlations 
with the other fit parameters. The strongest of these correlations are with ( and with R CT , 
which are illustrated in Figure HI 



5.3. Constraining Spin-Orbit Alignment 

We now use Model AC and the results of our MCMC analysis to examine the relationship 
between the spin axis of the black hole (the same as that of the jet; see Section [Tj), and the 
orbital angular momentum vector. We assume that the inclination of the orbital plane % is 



Gaussian distributed: % = 74?7 ± 3?8 flOrosz et all 120111 ). 



Our constraints on the locations of both axes are illustrated in Figure EJ which was 
derived using one million Monte-Carlo draws to represent each axis. This figure shows 
how readily our results are able to falsify the alignment hypothesis, even though we lack a 
measurement of the position angle of the orbital plane. Specifically, (1) over 80% of the sky, 
we are able to rule out the possibility that the spin and orbital axes are aligned^; and (2) 
the probability by random chance that the inclination angles agree so closely (see Fig. \5§ is 
less than 10%. 

Because the continuum-fitting method depends only on inclination angle (and not posi- 
tion angle), and because the difference between the inclination 9 of the spin/jet axis and the 
inclination i of the orbital plane is of critical importance in measuring the spin of J1550, we 
now use Model AC to determine 9 — i. Our results are shown in Figure E] where it is obvious 
that there is no evidence for any misalignment along the line of sight. That is, our results 
are consistent with 9 = i. We place upper limits on the absolute difference between orbital 
and spin inclinations of 8° and 12° at the 68% and 90% levels of confidence, respectively. 

Given the < 10% a priori chance that the inclination angles agree as closely as measured, 
our results provide support for the hypothesis that the two axes are aligned. However, 
without knowledge of the position angle of the binary axis (which can lie anywhere along 
the grey band in Fig. [5]), we cannot conclude whether they are, in fact, aligned. 



at 90% confidence 
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Table 2. Kinematic Model Settings and Fit Results 



Parameter Range Prior Shape a Model SI Model S2 Model AC 



9 (degrees) 


- 


89.99 


F 


53.9 ±0.7 


58.16 ± 1.7 


70.8+H 


r 

E h (10 45 erg) 


1 - 

io- 10 


1000 
-10 10 


LF 
LF 


21 q+390 

91.8+Jj 


crn+320 
OU_ 43 


36+1° 
c q+3.6 


D (kpc) 


3 


- T 


N(4.38i°f) 


3.07 ±0.06 


4 sn+ - 29 


4 48+ ' 43 

^■^°-0.34 


R cr (pc) 





-5 


LF 




0.46 ±0.03 


0.63 ±0.06 


5 


0.1 - 


- 10 4 


LF 




9401S8 


1041™ 


c 


io- 2 


- 10 2 


LF (maxfCC 1 ]) 






0.78 ±0.03 


min (Xp/M 








67.93 (543.4/8) 


41.59 (249.6/6) 


1.44 (7.21/5) 


min(x 2 / v) 








67.61 


40.94 


1.22 



Note. — The values quoted are the median parameter and symmetric 68% confidence interval (lex 
equivalent) derived from the MCMC run (as opposed to the single best-fit values). 

a F is flat, LF is log-flat, and N is a normal distribution. 

b Assumes tt-ism = 1 cm -3 and = 1°. 



c The lower bound on distance is taken from lHannikainen et al.l (120091 ) and the upper bound is 
derived using D < ~^== ( Mirabel fc Rodriguez 19991 ). 
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6. Radio Intensities and Asymmetric- Jet Models 

We now consider the radio intensity measurements discussed in Section [2] in order (1) 
to identify any intrinsic asymmetry in the jets and (2) to check the consistency of our 
kinematic model. In doing this, we are motivated by observations of the microquasar GRO 
J 1655-40, which in 1994 displayed multiple ejection events, each of which expanded and 
decayed on a time scale of a few days. The approaching and receding jets were found to 
be intr insically asymmetric; add i tionally, the sense of the asy mmetry changed from event to 



event (jHjellming fc Rupenl Il995t iMirabel fc Rodriguez! Il999l ) . 



Before introducing intrinsic jet asymmetry into the model, we first proceed by extend- 
ing Model AC to create Model RAC. This latter model retains the traits of Model AC 
but now incorporates Eqn. [8] and uses the additional free parameter A to model the radio 
data. The data set for Mo del RAC is likewise exte nded and includes its two radio intensity 



measurements (Section [21 lHannikainen et al.l I2009T ). Re- fitting the data using Model RAC 



and comparing with the results obtained for Model AC, we find a slight (2°) increase in 
the jet angle and similar small changes in the other parameters (Table [3]). The fit is good, 
X^/z/ = 1.44, and A, the decay rate of the jet emission, is positive and in the range ~ 1 — 5. 

We now examine intrinsic jet asymmetry, and introduce Model RAJ, a model that 
considers both kinematics and radio emission. In this case, the cavity is presumed to be 
symmetric, while the energy term (Eo/nQ 2 ) is allowed to vary between the eastern and 
western jets. Specifically, we set £ = 1, free r], and introduce the parameter q, which 
characterizes the type of asymmetry in the jets (Eqn. [8]). The fit results for Models RAC 
and RAJ are shown in Table [31 Model RAJ returns a significantly larger jet inclination angle 
than Model RAC, 9 ~ 82°, and it implies a large difference between the eastern and western 
jets, i]^ 1 ~ 15. Because q ps 0, for this model the gross asymmetry can be attributed to an 
east-west difference in the gas density (rather than an asymmetry in the energies or opening 
angles of the jets; see Section [3]). 

To assess the performance of Model RAC relative to Model RAJ, we exploit the sim- 
ilarities in the way these models are structured. In particular, their respective priors have 
identical form. Therefore, because we attribute equal likelihood to either type of asym- 
metry, we can apply the penalty normalization from Model RAC to Model RAJ. This 
yields the goodness-of-fit results shown in Table [31 Model RAJ is effectively ruled out: 
min(x2 PiRA j)-min(x2p iR Ac) = 8.2. 

We now test the strength of this result by considering a kinematic-only variant of this 
asymmetric-jet model, Model AJ, which ignores the radio intensity data. Model AJ has 
the virtue that it can be directly compared with our primary model, Model AC, because 
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both models have the same number of parameters (seven) and their priors are identically 
structured. As in the comparison above, we apply the penalty normalization of Model AC 
(Section [5]) to Model AJ. The fit results for the two models are given respectively in Tables [2] 
and [31 Based on the substantial difference in x2 p , min(x2 Pi Aj)— mm (x2 p ,Ac) = 7.5, and the 
even larger difference obtained when the radio-intensity data are included, we conclude that 
the asymmetric cavity model is favored over the asymmetric jet model at the 99% level of 
confidence. 

Unlike the manifestly asymmetric jets of GRO J1655-40, the available evidence indicates 
that the jets of J1550 are likely intrinsically symmetric: Model AC is favored over Model 
AJ, and Model RAJ implies an implausibly large (factor of 15) difference in the density of 
the ISM from west to east. 

In comparison with Model AJ or RAJ, our adopted Model AC gives a reasonable and 
satisfying description of J1550 as a system comprised of intrinsically symmetric jets prop- 
agating through an evacuated cavity with eastern and western walls located out at 0.6 pc 
and 0.5 pc, respectively. 



7. Discussion 



If we assume t hat the jets were pr oduced continuously over the day-long Eddington- 



limited X-ray flare ( jSteiner et all 120111 ). then the nominal total jet energy of ~ 10 46 erg 



implies that a significant fraction of the mass accreted onto J1550 during the flare was 
directly used to fuel the jets. Roughly, the initial mass in the jets was then ~ 10 24 g and the 
matter was accelerated to T ~ 10. 

We note that the moderate asymmetry we find (with the western cavity ~ 20% smaller in 



radius than the eastern) is opposite in sense from the asymmetry determined by lHao fc Zhang 
( 120091 ). We attribute this difference to several factors: Hao & Zhang simply adopted reason- 
able, ad-hoc values for several key parameters (8, E, and r ), and they found a high degree 
of asymmetry with t]^ 1 30 and ( = 1.4. (We note that this particular pair of values of 77 
and ( allowed a reasonable fit to be achieved to their data set.) By improving the quality 
and quantity of the astrometric data, we were able to determine that just one asymmetry 
parameter is required to explain the data, and that the resultant asymmetry is less extreme. 

Based on r esults obtaine d for the sub-pc scale (< 0.1 pc) jets of GRS 1915+105 and 
GRO J1655-40, iHeinzl (120021 ) has proposed that black hole microquasars preferentially in- 
habit environments that are under-dense compared to their AGN counterparts. Heinz offers 
several explanations, notably that microquasars may produce self-encasing low density bub- 
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Table 3. Additional Model Results 



Parameter 


Model AJ 


Model RAC 


Model RAJ 


9 (degrees) 


oc 0+2.4 


79 Q+7.4 


Q1 Q+5.1 


r 

E a (10 45 erg) 
D (kpc) 
R CI (pc) 
5 

C 
r/ b 

A c 
q d 


99+270 
zz -19 
91 q+83 
zio -65 

4.83 ±0.36 
0.46 ±0.05 

c-i n+1700 
Oiu -410 

0.065 ±0.014 


3711° 

6.1±|;1 
4 40+0.43 

^•^ 3 -0.35 

0.63 ±0.06 

yo -30 

0.78 ±0.03 
1.9±f;? 


i-4i±8:S 
soiSS 

q c;7+0.50 
00 1 -0.44 

S5+ - 04 

7/ln +3300 
'+ u -590 

0.068±g-~| 
-0.28l°;| 


min (Xp/^) 
min(x 2 /z/) 


2.95 (14.74/5) c 
1.11 


1.44 (8.63/6) 
1.11 


3.36 (16.81/5) c 
1.31 



Note. — The values quoted are the median parameter and sym- 
metric 68% confidence interval (ler equivalent) derived from the 
MCMC run. 

a Assumes nisM = 1 cm -3 and = 1°. 

b The form of the prior for £ and r] are identical (see Table |2j). 

C A flat prior is used for both A and q. The former is allowed 
to take values between [-10,10] and the latter is constrained to the 
range [-1,1.5]. 

d The penalty normalization for Model AJ is taken from Model 
AC; likewise that for Model RAJ is from Model RAC. 
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bles either as a remnant of the birthing supernova explosion, or via persistent kinetic outflows 
from the compact source. 

The enthalpy of the low density cavity in J1550, ~ 10 40 — 10 42 erg, is likely maintained 
by the steady (or quasi-stea dy) AU-scale jets known to be present in the har d or quiescent 
state of black hole binaries jRemillard fc McClintockl hood : iGallo et al.ll2006h . The ~ 20% 
measured asymmetry in the east-west extent of the cavity is unlikely to be a result of a high 
proper motion of the binary because this would require an extreme velocity ~ 0.1c. Rather, 
this asymmetry is easily explained as arising from a moderate 20% variation in the density of 
the ISM across the pc-scale region spanned by the jets. This supposition is quite plausible, 
given that J1550 is located only ~ 140 pc from the Galactic plane. 

One interesting feature of our best-fitting model is shown in Figure [21 The onset of 
X-ray emission for the western jet is first observed after the jet has reached the outer wall of 
the cavity, whereas for the eastern jet it occurs well before reaching the outer wall. Although 
there is not enough data to draw a firm conclusion, this difference in behavior suggests that 
our model oversimplifies by describing a succession of low-grade density jumps (from previous 
episodes of jet activity) as one single jump at R CT . Alternatively, perhaps one or several dense 
filaments of gas breached the eastern cavity walls, causing X-ray brightening at the shock 
front, but without contributing appreciable mass. 

We close our discussion by noting again that our lack of knowledge of the position angle 
of the binary restricts us to testing for spin-orbit alignment along the line of sight. The 
test we have performed nevertheless provides important support for the continuum-fitting 
measurement of J1550's spin, which u sed the orbital inc lination angle as a proxy for the 
inclination of the black- hole spin axis ( jSteiner et al.ll201ll ). For the case of J1550, we have 
shown that these two inclination angles are consistent within several degrees. 



8. Conclusions 



Building on earlier work by lHao fc Zhang! (120091 ) and IWang et al.l (120031 ) , we have used 
Chandra and radio imaging data to model the ballistic motion of the jets of J1550. We take 
the time of J1550's giant X-ray flare, which was promptly accompanied by the ejection of 
small-scale (~ 1000 AU) relativistic radio jets, as the launch date of the large-scale ballistic 
jets. Using our MCMC code and a kinematic model of the jets, we find that J1550 is 
enclosed in a pc-scale cavity that is moderately asymmetric, and that the jets are inclined 
by between 64° and 83° to our line of sight (90% confidence). These impulsive jets are 

extremely energetic, having been launched with a total energy of ~ 10 46 erg ( i deg ) ■ 
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By comparing our derived inclination angle for the spin axis of the black hole (taken to 
be the jet inclination angle) to the orbital inclination angle, we arrived at our primary result: 
We find no evidence for misalignment in our comparison of orbital and jet inclinations, and 
we conclude that the spin and orbital inclinations differ by < 12 degrees (90% confidence). 
This result has a likelihood of less than 10% of occurring by chance. 

Theory predicts that accretion torques acting over time will have brought most black 
holes into alignment with the orbital plane of their binary hosts. This prediction underpins 
the continuum-fitting method. In the case of J1550, our results provide support for such 
alignment and for the measured spin of its black hole primary. 
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40 F d 




2500 



Days since flare (MJD - 51076) 



Fig. 2. — The best-fitting model and fit residuals for the eastern jet (filled circles) and 
western jet (open circles). The cavity locations are marked by dashed horizontal lines, which 
indicate that the western wall (for the receding jet) is closer to the black hole than the eastern 
wall. For clarity, residuals for the coincidentally detected eastern and western jets are shown 
slightly offset in time. In the top panel, the error bars are smaller than the symbols. 
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Fig. 3. — MCMC results for Model AC. Probability densities are shown for each parameter 
on an arbitrary scale and have been obtained by marginalizing over all other parameters. An 
overlay for each prior shape is drawn as a dashed line. Note that the only two parameters 
which closely track the prior function are the system distance and Tq (at high values only). 
Otherwise, the prior contributes minimally to the parameter distribution. 
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Fig. 4. — Shown for Model AC are the MCMC density contours for the two parameters that 
correlate most strongly with inclination: cavity size R CI and the cavity asymmetry parameter 
(. The densities are calculated by marginalizing over all unshown parameters. Red contours 
mark the 68%, 95%, and 99.7% confidence regions about the most likely value, which is 
normalized to a density of unity. The central value of 9 changes from 70° to 80° as R cr varies 
from 0.65 pc to 0.58 pc, and as ( decreases from 0.78 to 0.74. 
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Observer's Line of Sight 
A 





Fig. 5. — The celestial sphere centered about J1550 with the observer situated along the 
pole. The angular momentum of the orbital plane is constrained to lie along the grey band 
(drawn with la width), and the spin angular momentum axis derived from the jets is overlaid 
in red (la about the most likely value). The position angle is completely unbounded for the 
orbital angular momentum, whereas the jets provide a tight constraint on the position angle 
of the black hole's spin axis (Section EJ). In fact, the uncertainty in the position angle of the 
jets is so small (±0.3 deg) that for purposes of illustration it has been tripled to make it 
visible in this figure. 

3000 E ; d 




I9-/I (degrees) 

Fig. 6. — The difference in orbital and jet inclination angles derived from the MCMC run 
of Model AC. The results show no sign of a misalignment along the line of sight; 68% and 
90% upper limits on the difference between inclinations are 8° and 12°, respectively. 



